Thermal Noise on Adiabatic Quantum Computation 
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The success of adiabatic quantum computation (AQC) depends crucially on the ability to maintain 
the quantum computer in the ground state of the evolution Hamiltonian. The computation process 
has to be sufficiently slow as restricted by the minimal energy gap. However, at finite temperatures, 
it might need to be fast enough to avoid thermal excitations. The question is, how fast does it need 
to be? The structure of evolution Hamiltonians for AQC is generally too complicated for answering 
this question. Here we model an adiabatic quantum computer as a (parametrically driven) harmonic 
oscillator. The advantages of this model are (1) it offers high ffexibility for quantitative analysis on 
the thermal effect, (2) the results qualitatively agree with previous numerical calculation, and (3) it 
could be experimentally verified with quantum electronic circuits. 

PACS numbers: 03.65.Yz, 03.67.Lx 
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Adiabatic quantum computation (AQC) offers an al- 
ternative route for achieving computational goals [1], 
compared with the "standard model" of quantum compu- 
tation based on the gate model. The basic idea of AQC 
is very simple: to maintain the system (computer) to 
stay at the (assumed unique) ground state with respect 
to a time-dependent Hamiltonian. For an isolated quan- 
tum system, this is in principle guaranteed by the quan- 
tum adiabatic theorem for sufficiently slow evolution, de- 
pending on the energy gap between the (instantaneous) 
ground state and the first excited state. Practically, AQC 
would be operated under some finite temperature which 
may not necessarily be negligible (recall energy gaps usu- 
ally shrink with the increase of the problem size) com- 
pared with the energy gap. The effect of thermal noise 
would then play an important role in determining the 
performance of AQC. Physically, the relaxation process 
(i.e. excitation to higher energy states) takes finite time 
to complete. Hence, we expect that AQC needs to be 
sufficiently fast as well. This sets another time scale due 
to the environment. Consequently, unless thermalization 
is not an issue, AQC would work only if the computation 
time lies within these two time scales. The question is, 
how to determine the latter time scale involved? We shall 
answer this question in this letter. 

This work is motivated by recent studies [2, 3, 4, 5, 
6, 7, 8, 9] related to the question of robustness of AQC. 
Concerning the noise effect on AQC, some of the models 
are based on either qualitative or perturbative arguments 
which are not verified by independent numerical investi- 
gation. Some of them are formulated in terms of parame- 
ters which are inaccessible experimentally. On the other 
hand, it was believed [7, 8, 9] that two-level approxi- 
mation would be valid for AQC, even if large number of 
excited states would be involved when the minimal gap is 
smaller than the temperature. It is therefore still unclear 
how "robust" AQC is against thermal noise. 

With these problems in mind, our goal here is to study 
the thermalization problem of AQC by modeling a Har- 



monic oscillator as a quantum computer. This model not 
only provides us with enough (infinite) excited states but 
also allows quantitative analysis. As we shall see, it could 
not be modeled by the two-level approximation. More- 
over, we shall quantify the effects on the performance 
of AQC through physical quantities such as temperature 
T, relaxation time I/7, energy gap A and computation 
time r* . The "anomaly" of this model may seem to be the 
evenly distributed energy levels. To verify the validity, 
we have compared it with the numerical simulation by 
Childs et al. [2], and found that the predictions of this 
harmonic model qualitatively agrees with their results. 
Lastly, this model is testable with the current quantum 
electronic technologies, e.g. simple RLC circuits. 

Adiabatic Quantum Computation — To define our adi- 
abatic quantum computer, there are only two adjustable 
and time- varying parameters, namely the "mass" 7714 = 
m [t) and the "spring constant" kt = k (t) . The time de- 
pendence of these two parameters, at this stage, is com- 
pletely arbitrary and is designed to sim,ulate (e.g. see ex- 
ample HI below) an adiabatic quantum computer. The 
(computational) system Hamiltonian Hs(t) is described 
by that of a standard parametrically driven harmonic os- 
cillator: 



Hs (t) = 



P 
2mt 



1. 
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which is associated with a set of (instantaneous) en- 
ergy eigenstates \nt) , with n = 0,1,2,..., satisfying 
the eigenvalue equation: Hs (i) \nt) — E^ (i) \nt). Here 
En (t) ~ {n + 1/2) A (t) is the (instantaneous) eigen- 
energy for the state \nt). The energy gap A (t) = 
yjktfmt ~ En+i (t) — En (i) does not depend on n, by 
definition. The initial state is assumed to be the ground 
state |0t=o) of Hs {t = 0). In the absence of the heat 
bath, the final state is given by U {t = tf) lOt^o); where 
U{t) = Texp{-iJ^Hs{t')dt') (with a = 1) is a time- 
ordered series. 

The computation is considered to be fail if the final 



state deviates significantly (due to excitation to higher 
energy states) from the desired ground state \Ot=f) of 
Hs {t = tf). This is best quantified by the fidehty F = 
|(0t=/| [/(i/) |Ot=o}| • Here since our goal is to study 
the thermal effect from the environment, we assume that 
AQC in the absence of the heat bath can be achieved 
(almost) perfectly, i.e., i^ « 1; violation of this condition 
may be considered as perturbation. 

Under this condition (and to zeroth order in A(t)), 
we may write U (t) |no) = exp(— i J^^ En (i') dt') \nt), and 
hence a relation which is needed later: 



U\t)atU{t) 



exp 



i / A {t') dt' ao 



(2) 



where at = y^mtAt/2 {x + ip/nitAt) is the (instanta- 
neous) annihilation operator for Hs{t). 

Ground State Occupation — In the presence of a heat 
bath, a mixed state representation p{t), or the reduced 
density matrix ps (t) = Trs {p {t)}^ is needed. The per- 
formance of the quantum computer is determined by the 
ground state occupation Pg = (Ot| ps (t) \0t), and in the 
coordinate space 



Pn 



dx'dx {x\ Ps (t) \x') ip* (x) (fit {x') 



(3) 



where ipt{x) = {x\ Ot) = (mtUJt) ex.p(^—mtUJtX^/2) is 
the (instantaneous) ground state wavefunction of Hs (t) . 
Before going into the technical details of the calcula- 
tions for PS (t), we first argue that, subject to the con- 
straints (a), (b) and (c) described below, the relevant 
quantity here is only the physical observable {x (t) ) = 
Tr{x^p{t)} (or the current fluctuation (/(i)^) in RLC 
circuits). 

The imposed constraints are (a) the heat bath can 
be approximated by a set of harmonic oscillators Hb = 
J2k ^^kbkbk, (b) the system-bath coupling Hsb is bilin- 
ear e.g. terms like x{bk + f'J,), and (c) the initial state of 
the bath is in a thermal state ps = e^'^-^^/Tr [e~'^"''] 
and the system is in the ground state of Hs (0), i.e., 
p{0) = |0t=o) (Ot=o| ® Pb- To proceed, we write 



\x') {x\ 



1 r°° 



,i{fJ,p+i^x) 



(4) 



where fj, = x — x' (and i^ is a just dummy variable). 
This form suggests that wc have to evaluate the quan- 
tity (e*(^P+''*') = rr{e*(w+''*)p(i)}, which is equal to 
exp(— ((/ip + i'x)^)/2) from the Bloch identity. Now, as 
verifiable by the master equation, here we claim that 
(flj) = (aj ) = 0. By completing the Gaussian integrals 
in Eq. (3) and (4), we finally arrive at a very compact 
form for P„: 



Pn 



1 



l + n{t) 



(5) 



where n{t) = {a\at) = Tr{a\atp{t)}. Thus, as adver- 
tised, (x^) = {h/2mtAt) {2{alat) + 1) is the only quan- 
tity needed to determine Pg. 

Master Equation — We shall obtain n{t) through the 
master equation approach [10]. Here the full Hamiltonian 
is divided into three parts: H = Hs {t)+HB+HsB where 
the first two terms have been defined. We assume that 
the coupling term Hsb is a time independent operator 
(i.e. indepdendent of the mass mt and spring constant 
kt of the oscillator), and is explicitly given by 



Hsb — 



k 



(6) 



Note that there could be a frequency renormalization 
(lamb shift type), which modifies the ground state wave- 
function. This effect would be small for weak damping 
A(i) > 7(t) where 7 (t) = ri{t)/m{t) (cf. Eq.(7) and 
(8)). Second, even for ohmic damping (time indepen- 
dent r] (t) = rj), the relaxation rate 7 (t) oc l/m (t) (or 
time to reach equilibrium) depends on the system pa- 
rameter (here the "inertia" m(t)), and therefore may be 
time-dependent . 

To continue, we shall keep the standard assumptions 
for the master equation, namely (i) product initial state, 
(ii) Born-Markov approximation (i.e. weak coupling and 
short memory time), and (iii) rotating wave approxima- 
tion (i.e. ignore fast oscillations). Subject to these con- 
straints, the master equation is given [10] by 



di^'^l^JO 



dt'TrB{[HsB{t), [HsB{t'),Ps{t) ® Pb]]}, 



where in the interaction picture: ps{t) = W{t)ps{t)U{t) 
and HsB{t) = U'' {t)HsBU{t). If we write x = 
y^h/2mtAt{at + aj), and from Eq. (2), we obtain in- 
teraction terms similar to that of an ordinary (i.e. with 
mass and spring constant fixed) harmonic oscillator, ex- 
cept the replacement: (a) Agt -^ J^ dt' A (t') and (b) 
toqAo — > mtAt- Consequently, the resonating modes 
LUk ~ A (t) would be time-dependent, and hence the fric- 
tion "coefficient" 77 (t) = 7rJ(At)/At, where J (oj) = 
Sfe5fc^('^fe — w), would also be a function of time i.e., 
with a classical equation of motion (neglect the frequency 
renormalization) : 

m{t)-^{x)+rj{t)-^^{x)+k{t){x)=0 . (7) 

The exception is the ohmic case, where J (lu) oc uj and 
hence 77 (t) = 770 is independent of the variation in the 
mass and the spring constant (e.g. RLC circuit). Finally, 
the equation for n (t) = Trjajaops (t)} is obtained from 
the master equation: 



dt 



n{t)^~^it)init)-Nit)) 



(8) 



l) . This is the key result 



where N {t) = 1/ (e^(*)/'=«'^ 
of this paper, since the performance of AQC is deter- 
mined entirely by n{t). Note that even for the case of 
ohmic damping, the relaxation rate 7 (i) =77 (t) /m (t) 
is in general time dependent. With the initial condition 
n (0) = 0, this equation can be solved numerically to 
obtain the ground state occupation Pg at time t. 

Although in our model the time dependence for the 
energy gap is completely arbitrary, for the purpose of 
understanding the structure of the thermal excitation we 
assume that the gap has a Landau-Zcncr type variation: 



A(t) 



A2 

max 



{l^t/r^r + Al,^ 



(9) 



Ainax, and 



where for A^ax > A,„i„, A (0) = A (2r,) 
A (t») = Amiii. Except near the region i w tq, the rate 
of change of the energy gap is Vs = Amax/T* . For sim- 
plicity, we shall consider the ohmic case only and assume 
that m (t) = mo is time- independent, which makes 7 
time-independent as well. 

The following examples are chosen to demonstrate re- 
spectively that: (I) when thermalization is important 
(i.e. A (i) < kgT), the computation speed Vs needs 
to be fast, compared with the "natural" speed of the 
bath Vb = ikgT. Wc quantify this by defining R = 
Vb/Vs = 7fc_B7'r,/Amax- (II) After passing through the 
gap minimum, when the energy gap is larger than the 
temperature, i.e. A (i > r*) > ksT, relaxation towards 
the ground state (increasing Pg) has the simple e"''* 
dependence, and in contrast with that in Rcf. [9], docs 
not depend on A^j^jj^/Amax ui the exponent. (Ill) This 
toy model qualitatively agrees with numerical calculation 
based on more realistic Hamiltonian. 



Example I 



Consider the case where A^ax ^ kgT 



and Anjin «C kgT. It is possible to approximate N (t) 
ksT/A {t). Substitute this into Eq. (8), and with 7 (t) 
7, we have 



nit) 



■yksT 



-ft 



ds 



{l-.s/T,y+e^ 



(10) 



where e = Amin/Amax *C 1. For 7r* < 1, the integrand 



is dominated near s ^ t^,. Taking e'*'''* 
grating explicitly, we have 

n(t) = A(t)i?e-'^(*-^*) 



and inte- 



(11) 



where X{t) = In 2 - ln[^Je^ + (1 - t/r^f + (1 - t/r^)], 
and R = jkBTr^ /Amax as defined above. Figure 1 shows 
that this expression for n(t) is in good agreement with 
the result by direct numerical integration for n(i). From 
Eq. (5), we conclude that the thermal effect is not impor- 
tant even if A (t) < ksT provided that 7t« <C 1. More 
precisely, we require i? <C 1, or V5 » Vb. This minimal 
speed limit for AQC could not be seen by the two-level 
approximation [9]. 
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FIG. 1: Simulation of AQC with harmonic oscillator under 
thermal noise. The a::-axis is rescaled time 'yt. 1(a) and 2(a), 
in unit of fcsT, show the energy profiles Eq. (9) for two classes 
of AQC (but with the same R). 1(b) and 2(b) show the dy- 
namics of the mean excitation of the oscillator. The solid lines 
are numerical integration from Eq. (8), and the dashed lines 
are analytic approximations Eq. (10) and Eq. (12). 1(c) and 
2(c) are the corresponding ground state probabilities obtained 
from 1(b) and 2(b), with the relation Eq. (5). For compar- 
ison, the dotted lines are the thermal equilibrium values of 
P - 1 _ „-/3^ 



Example II — Here we consider the possibility of re- 
laxation after passing through the minimum. In other 
words, we consider nit) when A (t > r*) > kgT, while 
Aniin < ksT . This situation should not be very common 
for AQC, as it suggests that thermalization from the heat 
bath would yield better performance. To start, wc could 
not invoke the same approximation as in example I. How- 
ever, as long as i? < 1, A^ (t) is still sharply peaked at 
i « T, . Base on this observation, skipping the details, we 
obtain an approximate solution which is valid only for 



;(t>n) = Ki?e-^(*-^-) 



(12) 



where k = 2a + 21n [fcsT/Amin (1 — -R)] and a = 
^^ e^^\n[2x)dx = 0.116. Figure 1 shows that this 
expression qualitatively agrees with the direct numer- 
ical integration for n{t) for t > T.^,. Again, we con- 
clude that R plays an important role for determining the 
performance of AQC. It is suggested [9] that the com- 
bination A^jijj/Amax would be important for the relax- 
ation process in the two-level approximation. We have 
tested with different ratios of Amin/Amax (while keeping 
Amin ^ A„iax), but WC did not find explicit dependence 



of Amin in the exponent of n{t). 

Example III — So far we have compared our results 
with that from the two-level approximation. Would a 
realistic Hamiltonian (with non-uniform distribution of 
energy gaps) for AQC, in some sense, look like a har- 
monic oscillator (uniform gaps)? If yes, then based on 
the results above, it may be possible to approximate the 
final ground state occupation by the formula: 



P. 



1 



l + aR 



(13) 



where R = jhsTT,:/ (Amax ~ Amin) is generalized to in- 
clude cases where Amin/Amax is not negligible, and a 
is a fitting parameter. For example 1, at t = 2r*, a = 
2e~''''^* In (2Aniax/Aniin), and for example II, a — Ke~^'^' . 
The former does not depend on T, and the latter depends 
on T weakly (logarithmically). 

In Ref. [2], an algorithm solving the so-called "three- 
bit exact cover" (ECS) problem, in which the energy gap 
(here taken as A (f)) between the ground and first excited 
state of Hs (t) is significantly larger than the rest when 
A (t) = Aniin- We extract, from FIG 2 of that paper, the 
final probability Pg and R, and estimate the correspond- 
ing a by the relation in Eq. (13), as suggested by our 
harmonic oscillator model. The results are shown in Ta- 
ble I. We found that both sets of data are consistent with 
the conjecture that Pg decreases with R. For data I, the 
fluctuation for the value of a is relatively small (about 
20% from the mean) , and the data point for high temper- 
ature case (fcsT/Aniax = 10) deviates significantly with 
the rest. This is anticipated from our experience in ex- 
amples I and II. For data II, the fluctuation is relatively 
larger (about 40% from the mean). This may due to that 
the ratio Amin/Amax — 0.425 is a bit too large for our 
simple formula Eq. (13) to be accurate. We conclude 
that this harmonic oscillator model may provide reason- 
able estimation for some realistic AQC problems. 

For AQC involving many more degrees of freedom, it 
may (either computationally or experimentally) be chal- 
lenging to obtain the eigenenergy spectrum. However, 
this model can still be applicable. We first determine 
the a and R for an AQC with a relatively small prob- 
lem size n (i.e. what was done in Table I). We then 
gradually increase n and determine the (average) scaling 
r„ < 1 of the energy spectrum (typically the first few 
lowest energy states are enough). Then, to estimate the 
same AQC (under the same temperature) with large N, 
the harmonic oscillator model suggests the ground state 
probability be given by the formula Eq. (13) with the 
replacement A^ax -^ Tw A,„ax (assuming r* and 7 are 
fixed). 

Conclusions — We have introduced a harmonic oscilla- 
tor model for a quantitative study on the effect of thermal 
noise on AQC. For ohmic damping, we showed that AQC 
is considered fast, if the combination R = "fkBTr^/A^i^y^ 
is small. This model suggests a simple relation for es- 



TABLE I: Simulation of AQC with Harmonic Oscillator. The 
data (Pg and R), taken from the numerical simulation (FIG 
2) of Childs ei al. [2], are fitted (a being the fitting parame- 
ter) with the formula in Eq. (13) suggested by the harmonic 
oscillator model. The standard deviation (excluding the last 
data point ksT = 10) of a for data I (II) is about 20% (40%) 
from the mean value 1.68 (0.81). 
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1/10 
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Data I: 


P, 


0.79 


0.53 


0.30 


0.15 


0.08 


A„,i„=0.301 


R 


0.14 


0.72 


1.43 


2.86 


14.3 




a 


1.86 


1.24 


1.63 


1.98 


0.80 


Data II: 


P, 


0.89 


0.70 


0.42 


0.19 


0.08 


A^i„=0.425 


R 


0.17 


0.87 


1.74 


3.48 


17.4 




a 


0.71 


0.49 


0.79 


1.23 


0.66 



kgT and A^in are in unit of Amax- 

timating the fidelity (cf. Eq.(13)) of general AQC; this 
relation qualitatively agrees with the previous numerical 
simulation. This model can also be verified with quan- 
tum RLC circuits, and therefore can act as a test bed 
for future theoretical and experimental investigation on 
AQC. 
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